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Relativistic electron transport through an oscillating barrier: 
wave packet generation and Fano-type resonances 
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Transport properties of massive Dirac particles are investigated through an oscillating barrier. 

The Floquet quasienergies related to the time-dependent potential appear both in transmission and 
reflection as sidebands around the incoming electron’s energy. We take all relevant sidebands into 
account and present time averaged transmission and reflection probabilities in a wide energy range. 
Most qualitative features of scattering on a static barrier - like Klein paradox - are still visible, 
but the transmission probability in the evanescent regime observably increases due to the oscillation 
of the potential. The strongly inelastic scattering process is shown to lead to multiple Fano-type 
resonances and temporal trapping of the particles inside the oscillating potential. We also present 
a detailed study of the time evolution of the wave packets generated in the scattering process. Our 
results can be relevant for graphene with an induced energy gap. 

PACS numbers: 73.23-b, 03.65.Pm, 72.80.Vp 


I. INTRODUCTION 

Quantum scattering by time-harmonic potentials is an 
important and vivid research area. It provides deep un¬ 
derstanding of a rich variety of interesting and partly un¬ 
usual phenomena in strongly driven quantum systems. 
Photon-assisted tunnelingii^ is a remarkable example 
showing that the presence of an alternating field can lead 
to strongly inelastic processes. Nonrelativistic quantum 
mechanical scattering on barriers with oscillating height 
has been investigated intensively - mainly in the con¬ 
text of traversal time and photon assisted transport, see 
e.g. Refs. It has also been proven that Floquet 

theory^ provides an efficient tool for the investigation 
of various time-dependent scattering processesi^^ Fano- 
type resonance a^d^ can appear in this context due to 
transitions between sideband states and bound states. 
Developments in the experimental techniques during the 
last decade enable that these results have the potential 
of direct applications in the rapidly expanding field of 
meso- and nanoscale quantum devices 

Much of the theoretical works published so far studied 
transport through time-dependent potentials in a nonrel¬ 
ativistic frameworkdr— Transport related prob¬ 
lems with oscillating spin-orbit interaction have been 
studied e.g. in Refs. |15l - [l7f . A few recent papers^^— 
treating massless Dirac particle scattering on time- 
harmonic potentials, are inspired by the unique elec¬ 
tronic dispersion relation of graphene . In this single 
layer of hexagonal carbon atoms - as it has been demon¬ 
strated experimentallyddr— - the carriers exhibit strik¬ 
ing relativistic features like Zitterbewegung;ddi21 Klein 
paradojjd^— and Klein tunneling.— idi 

Substrate-induced bandgap in epitaxially grown 
graphenedd opened the way for its usage as an electronic 
material. This induced bandgap leads to a finite mass for 
its charge carriers which obey massive Dirac equation, 
and the energy dispersion relation is no longer linear in 
momentum. Based on this, Klein tunneling of massive 


Dirac fermions through a static barrier^ and massive 
Dirac electron tunneling through a time-periodic poten¬ 
tial in single layer graphene^ were studied. 

In the present paper we investigate theoretically the 
relativistic scattering of massive Dirac particles on a 
time-periodic rectangular potential barrier (see FigllJ, 
using Floquet’s technique. The method we use is tech¬ 
nically similar to that of Ref. [s^, but we concentrate 
on different physical aspects of the problem by present¬ 
ing spacetime resolved results and discussing Fano-type 
resonances. Besides its relevance to graphene with a 
bandgap, this work may also contribute to the manip¬ 
ulation of relativistic charged particle beams by powerful 
laser pulses— “— 

Our paper is organized as follows: in Section II we 
define the model and derive the equations for the ampli¬ 
tudes, to be solved numerically, from the one-dimensional 
Dirac equation using Floquet theory. In Section III we 
present the physical contents of the numerical solution 
(e.g. Klein tunneling) in terms of the cycle-averaged and 
the time-dependent transmission and reflection probabil¬ 
ities. We explore the details of the scattering with the 
help of space and time dependent charge current and elec¬ 
tron density. Finally, we analyze and explain in detail the 
Fano-type resonances^iiS found in the transmission and 
reflection curves. We close our paper by summarizing the 
results and drawing the conclusions in Section IV. 


II. MODEL 

In this work we consider a two band model of graphene, 
where the bands are separated by a band gap. This is mo¬ 
tivated by the experiments?^ where the electronic struc¬ 
ture of this unique material has been modified and the 
degeneracy of the Dirac point at the intersection of the 
valence and the conduction bands could be removed by 
growing the sample on a SiC layer. As shown by the mea¬ 
surements, in this case the dispersion relation became 
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A. Solution of the Dirac equation with oscillating 
potential 

Inside the region 0 < z < L the Dirac equation reads 
= H{t)ip{z,t), (3) 

with 

H{t) = Ho-\-Vi{t) = ca-i +/3mc^ + Vo +Vi(i), 

(4) 

where the standard 03 and /3 matrices appear. Since 
[H{t),H{t')] = 0, an eigenstate 

Hoip = £ip (5) 


FIG. 1: Schematic view of the one dimensional scattering 
problem we consider. The harmonic oscillation of the poten¬ 
tial has an amplitude oi hfl and angular frequency of v. That 
is, V{t) = Vo h^l cos ut in region 2, while the potential is 
zero in regions 1 and 3. A monoenergetic electron wave is 
assumed to impinge the oscillating barrier inducing reflected 
and transmitted waves in region 1 and 3, respectively. The 
input energy Eq is always larger than mc^ (half of the induced 
bandgap), but the magnitudes of Eq, Vo and ftfl relative to 
each other were varied in our calculations. 


similar to that of a massive relativistic Dirac particle, 
with a band gap A (which appears as 2mc^ in the Dirac 
equation) as large as 0.26 eV. We consider a one dimen¬ 
sional model where this gap is influenced by a constant 
plus a harmonically varying potential: 


V(t) = Vq-\-V i{t) = Vq-\- h£l cos vt 


( 1 ) 


which is similar to the gapless model of Ref. . We also 
note that the nonrelativistic version of a similar problem 
has been considered in Refs. [ug in the context of mul¬ 
tiphoton ionization. Outside the oscillation region the 
potential is zero, see Fig.[TJ A monoenergetic spinpolar- 
ized free electron wave is assumed to impinge the oscillat¬ 
ing potential barrier, i.e, in the standard representation^ 
we have 
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where the direction of propagation has been chosen to be 
the z-axis. The spinor above is a solution of the Dirac 
equation with Hko = According to the 

geometry shown in Fig. [TJ we choose the positive sign 
here. 


can be used to construct a solution to Eq. ©: 




<^( 0 ) 


W)dr 


( 6 ) 


For a fixed value of k, (i.e., spatial dependence of 
we have 


£^{k) = ±\/-\- h?k‘^c^ Vo, (7) 

which are both doubly degenerate (due to the two pos¬ 
sible spin directions). Since the interaction is indepen¬ 
dent of spin [the terms Vb -I- Vi{t) are proportional to 
the unit matrix in Eq. ©], the solutions of Eq. © that 
correspond to the monoenergetic incoming spinor ([ 2 ]) as 
a boundary condition, have nonzero components at the 
same positions as Therefore it is sufficient to con¬ 
sider only 

^±(z,t) = e*'="u±(fc) (8) 
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( 10 ) 


that clearly satisfy Eq. ©. Using Eq. 0 , it is read¬ 
ily seen that the spinors above do not depend on Vq- 
Thus e.g. Eq. m we can be rewritten as ipin{z,t) = 
Eko^-i^t^+l^ko). [In this case £+{ko) = Eq.] 

Note that the only restriction in Eqs. 0 ( 101 ) concern- 
ing k is that h?k^ > —m?c^ [ensuring that £^ are real, 
see Eq. 0 ]- This means that evanescent solutions with 
purely imaginary k are also allowed. 
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Using the Jacobi-Anger identity 

oo „ 

( 11 ) 

n— — oo 

where Jn denote Bessel functions of the first kind, we see 
that the frequencies appearing in the time evolution are 
given by £^[k)/h + ni', with integer n. Note that since 
the differential operator given by Eq. (|3]) is periodic in 
time (T = 2TrIv), Floquet theory^i^ can be applied. The 
states ([5|) are orthogonal in the spinor sense, they can be 
considered as elements of a time-dependent basis. Apart 
from the factors exp{—i£^{k)t/h), these solutions are pe¬ 
riodic, thus can be called the (nonequivalent) Flo¬ 

quet quasienergies. The term nonequivalent means here 
that e.g. £~ (fc) = £~ (fc) -I- nhi' can also play the role of a 
Floquet quasienergy (with n being an integer), but states 
with time-dependences exp[—i£~{k)t/h — i{£l/i')sm{i't)] 
and exp[—i£~{k)t/h — i(£l/ly) + inhiy] are dynam¬ 

ically equivalent. On the other hand, £'^{k) and £~{k) 
correspond to qualitatively different dynamical behavior 
(unless their difference is an integer multiple of hv). 


B. Fitting the solutions 

In the previous subsection plane wave solutions of the 
Dirac equation were obtained. It was shown that for a 
given (real or purely imaginary) value of fc, Eq. ([3|) is sat¬ 
isfied by the spinors given by Eq. (jS]). However, 

as we shall see, in order to obtain a solution to the prob¬ 
lem over the whole 2 -axis, several (in principle infinite 
number of) different wave vectors are needed. 

According to the previous subsection, in region 2 
(where the potential oscillates) whenever a frequency £/h 
appears in the time evolution, the harmonics £ jh nv 
are also present (n = ..., —2, —1,0,1,2,...). However, if 
we want to impose continuity of the spinor valued wave 
function at the boundaries 2 = 0 and 2 = L, the linear 
equations have nontrivial solutions only if the input fre¬ 
quency equals to one of the harmonics mentioned above. 
In other words, the frequencies we have to take into ac¬ 
count are 


uJn = En/h = Eo/h +niy, (12) 

where Eq is the well defined energy of the input spinor, 
see Eq. Q . In region 1, the only right propagating spinor 
(see Fig.[T]) is the input; a particular solution of the Dirac 
equation corresponding to frequencies is given by: 


LJti >0 

+ E (13) 

CDn <0 


where 
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j m^c^ — E^ 
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if < w?c^. 


(14) 


The signs here have been chosen such that the terms 
proportional to Vn in Eq. (fT^ describe either reflected or 
evanescent waves (with exponentially decaying amplitude 
as 2 —>■ — 00 .) Analogously, in region 3: 


a;n >0 

+ E e*'=""u-(fc„) (15) 

CJti<0 


The problem is somewhat more complicated in region 
2 , where 
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VqY > m?c^ 
VqY < rn^c^ 


(16) 


are the solutions of 


£^{kn)=En. (17) 

However, in this case, due to the oscillation of the poten¬ 
tial, a given wave number k'„ corresponds to not only a 
single ujn- For the sake of simplicity, we collect the coef¬ 
ficients of the two types of spinors given by Eqs. (igiTUl) 
separately, and write 

^'2(^,0 = ^'J(2:,<) + ^'^(2:,^)- (18) 

The first term here is given by 

X [a„ e<^u+(A:;) + 6„ e-<^u+(-fc;)] , 

(19) 

where, as we can see, both propagation directions appear. 
Additionally, due to the finite size of the region, exponen¬ 
tially growing spatial dependence is also allowed. Next 
we insert the Jacobi-Anger expansion (El) and obtain an 
equation where the frequencies appear explicitly. For 
the sake of brevity, we present this step only for the terms 
containing the spinor u~: 

'^2 {z,t) = E E ® 

X [a„ e<^n-(/c;)+6„ e-<^w-(-fc;)' . 

( 20 ) 

The condition of continuity at 2 = 0 (2 = L) now 
can be formulated by evaluating 4'i(0,t) and 4'2(0,t) 
[4'2(T,t) and 4'3(L,t)]. Working in frequency domain. 
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as an example, the contribution of ^'2 {z = 0) to the 
frequency component uji is given by: 

(“) + bnu~{-k'j] . ( 21 ) 

By comparing the coefficients of each w; at the bound¬ 
aries, we obtain an infinite system of linear equations for 
the unknown coefficients {rmtman,bn}. However, since 
the Bessel functions J for a given argument generally 
decrease as a function of their index, correct numeri¬ 
cal solutions could be obtained by taking only a finite 
number of frequencies into account. If we consider a set 
{ojn '■ n = —N ,..., 0,... fV}, there will be 8N + 4 fitting 
equations and the same number of unknowns. As we shall 
see in the next section, the time averages of the transmis¬ 
sion and reflection probabilities provide an efficient tool 
for monitoring the accuracy of the numerical method: if 
their sum is not as close to unity as we require, N has to 
be increased. 


III. RESULTS 


The time-dependent reflection and transmission prob¬ 
abilities are given by the ratio of the transmitted and 
reflected currents to the incoming one: 


m 


jrjt) 

jin 


T{t) = 


jtjt) 

jin 


Eq + mc^ 
2chkQ 
Eq + mc^ 
2chko 


4'Ja34'3(L,t), 


( 22 ) 


where 


= ^1 - V'in, (23) 

i.e., the reflected (or exponentially decaying) part of 'I'l. 
Let us note that the parameter ranges we use in the fol¬ 
lowing are ideal to see and identify the physical processes 
that are responsible for the effects to be presented. To 
this end we use ’’natural units” (i.e., h = l,m = 1 and 
c = 1) for the figures. 


A. Cycle-averaged reflection and transmission 
probabilities 

Since the time dependence of both T and R contains 
factors exp[—i(a;„—Wm)t] = exp[—w(n—m)t], these func¬ 
tions are periodic, T{t+T) = T{t), R{t+T) = R{t), with 
T = 27r/i/. First we consider the time average of the re¬ 
flection and transmission probabilities, 

(T) = f T{t)dt, (R) = f R{t)dt. (24) 
Jo Jo 



FIG. 2: (T) as a function of barrier height Vo for the indicated 
values of the oscillation amplitude Q. Additional parameters: 
Eo = = 0.2, L — 10. The physical meaning of the 

spikes will be discussed later in Sec. Ill C. 


Using Eqs. m cllld (|15l) , WG hcLVG 
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(25) 


Note that since the scattering problem is periodic in time, 
the cycle-average of the incoming current should be equal 
to (jt) - Or) = \{jt) \ + I Or) I • In other words. 


(T) + (i?) = 1 (26) 


should hold. This requirement can be used to monitor the 
accuracy of our calculations. Since we truncate the infi¬ 
nite system of equations, it is not necessary, that Eq. (ESD 
is satisfied. However, if the populations of the states that 
are neglected due to the truncation are negligible, the er¬ 
ror can be kept below an acceptable limit. Eor the results 
to be presented in the following, |1 — (T) — (i?)| < 10“®, 
and to achieve this limit it was sufficient to truncate the 
system at A = 20 — 25 depending on various parameters. 


Figs. [2]l4] show (T) as a function of Vq for the weakly 
relativistic, relativistic and ultrarelativistic cases (when 
Eq is close to mc^, Eq = 2mc^ and Aq = lOmc^, re¬ 
spectively). The dashed black curve - as a reference - 
corresponds to the case of H = 0 (non-oscillating barrier) 
in all figures. The most important point that Figs. [2]l4] 
have in common is that for U = 0 (and oscillations with 
small amplitude) (T) is practically unity when Vq <C Eg; 
it is almost zero when Eq — mc^ < Vq < Eq + mc^ 
and converges (in an oscillating way) to I again, when 
Vq > Eq + mc^. This is a well known behavior, that can 
be understood readily by investigating Eqs. (fTBl) : The 
solutions in region 2 (see Eig. [T]) are propagating waves 
in the first case, real exponentials in the evanescent do¬ 
main {Eq — mc^ < Vq < Eq + mc^) and propagating 
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FIG. 3: Panel (a): {T) as a function of barrier height Vo for 
the indicated values of the oscillation amplitude fl. Addi¬ 
tional parameters: Eq = 2m (?= 0.2,1/ = 10. Panel (b): 
Zoom into the interval Eo — rru? < Vb < So + me? of panel 
(a). 


FIG. 4: (T) as a function of barrier height Vo for the indicated 
values of the oscillation amplitude Q. Additional parameters 
in panel (a): Eq = = 0.2, L = 10, panel (b): Eq = 

lOmc^, u = 0.2, L = 1. 


waves again when Vo > Eq + mc^. In other words, (T) 
reproduces the Klein paradox^^Sr^ for small values of 17. 

The oscillations that can be seen in Figs. [2]l4] are sig¬ 
natures of quantum mechanical interference: there are 
maxima (minima) in (T) when the spinor valued waves 
interfere constructively (destructively) at z = L. When 
the ratio 17/;/ is small, terms corresponding to wavenum¬ 
bers ko and /cq dominate the dynamics. For larger val¬ 
ues of 17/;/, the expansion (fTTI) in terms of Bessel func¬ 
tions contains numerous frequency components resulting 
in more complex oscillation patterns in Figs. [2E1 

The figures show an additional, important effect, 
namely the gradual disappearance of the pronounced flat 
minimum of (T) (Vq) as either L is decreased, or 17 in¬ 
creased. The first case is related to the role of the 
evanescent solutions, since tunneling becomes increas¬ 
ingly efficient when the width of region 2 is decreased. 
When the amplitude of the potential oscillations is in¬ 
creased, more and more frequency components play a 
relevant role in the dynamics. Some of them corresponds 
to (quasi)energies En that are higher than the oscillating 
barrier, and consequently the related part of the spinor 
valued waves are transmitted with a high probability. 


This effect is still present in the case when the energy 
of the incoming spinor, Eq, is below the minimum of the 
oscillating potential: 


Eo-mc^ <Vo-hVt. (27) 


Without oscillations, (T) = T would be practically zero 
in this case. Fig. [5] corresponds to parameter values 
where Eq. (EID is satisfied and T < 10 ^ when 17 = 0. 
However, as we can see in the figure, orders of mag¬ 
nitude higher cycle-averaged transmission probabilities 
arise when 17 = 0.3 (in natural units.) The dependence 
of (T) on the frequency v shown in Fig. [5] tells us that 
although lower values of v with a fixed 17 means that 
more frequencies uJn should be taken into account, this 
effect is weaker than the fact that higher values of v cor¬ 
responds to larger steps in the ladder Wn = wq + nv- That 
is, as a tendency, (T) [v) is an increasing function in the 
parameter range given by Eq. (123. 
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FIG. 5: (T) as a function of i/, for parameters Eo = 

3mc^,Vo = 2.75, L = 10,17 = 0.3. Note that the transmis¬ 
sion probability in the static case for the same parameters 
(evanescent solutions) has the order of magnitude of 10“®. 
Note that the physical relevance of the abrupt changes seen 
in this figure will be analyzed in detail in subsection IIII Cl 


B. Spacetime dependence of the process: wave 
packet generation and propagation 

It is instructive to investigate the space and time de¬ 
pendent quantities 

p(z, t) = t), j{z, t) = c'^/\z, t)a 3 '^/{z, t) 

(28) 

that satisfy the continuity equation 

§-^piz,t) = -^^j{z,t). (29) 

The spinor vk in Eq. (1^ above stands for 'hi, f = 1, 2, 3, 
depending on whether z is in region 1, 2 or 3, respec¬ 
tively. Note that for = 0, p does not depend on time, 
and consequently j is constant as a function of z. Consid¬ 
ering an oscillating potential, the spacetime dependence 
of both quantities are considerably more interesting, and 
sheds light on the cycle-averaged results presented in the 
previous subsection. 

The physical phenomena being responsible for the re¬ 
sults presented so far are most visible by focusing on 
parameters yielding (T) = (R) = 1/2. Figs. H] and [7] 
show p{z,t) and j{z,t) in this case, for different values 
of n. As previously, the patterns that can be seen are 
more complex for larger amplitude of the potential oscil¬ 
lations, but the qualitative features of Figs. [5] and [7] are 
the same. Considering the electron density, it has certain 
pronounced maxima in region 2 (the value of p at these 
points can be larger than anywhere else by a factor of 
two), i.e., we observe a kind of temporal ’’trapping” of 
the population inside the oscillating barrier. 

We can also see structured wave packets that propa¬ 
gate in the positive (negative) z direction in region 3(1) 
due to the ’’pumping” (see e.g. Refs. caused 

by the oscillating barrier. Considering the time intervals 
when these wave packets are released, it is instructive to 



t v/27t 


FIG. 6: Density plot of p{z,t) in panel (a), and of j{z,t) in 
panel (b), and plot of Eq — me? — V{t) in panel (c), for the 
parameters L = 10, Vo = 1.95mc^, Eq = 3mc?, and D = 0.2. 
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FIG. 7: Density plot of p{z,t) in panel (a), and of j{z,t) in 
panel (b), and plot of Eq — me? — V{t) in panel (c), for the 
parameters L = 10, Vo = 1.85mc^, Eo = 3mc?, and D = 0.5. 


see panel c) in Figs. Eland [71 where Eo—mc^—V(t) is plot¬ 
ted as a function of time. When this quantity is positive, 
the potential barrier is lower than Eq, and we see wave 
packets leaving the central region in the forward direction 
(transmission). Oppositely, when Eq — mc^ — y{t) < 0, 
we have mainly reflection. 
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The current density’s spacetime behavior resembles 
that of p, but in this case the sign contains an addi¬ 
tional information. In region 3, j is by construction 
always positive, although its magnitude varies. In this 
region p and j have very similar spacetime dependence. 
For z < 0, however, both the incoming and the reflected 
spinors contribute to the current density, and their inter¬ 
ference can result in negative or positive j (depending on 
whether or is the dominant, respectively.) Note 
that when {R) Ki 1, j{z,t) is practically zero except re¬ 
gion 1, where it represents a truly alternating current, 
with max[j(t)] Ri — min[j(t)]. Obviously, in this case the 
cycle average of the current is zero everywhere. 

C. Fano-type resonances 

As it is known, whenever the energy of a scattering 
state coincides with that of a bound state, or in other 
words, these two different kinds of eigenstates happen 
to belong to the same degenerate energy level, then the 
transition probability between the states becomes large 
and a resonance occurs in the transmission. In general, 
this effect is known as a Fano resonance<2ii^ 

The role of Fano-type resonances has been discussed 
earlier in the context of the model of gapless graphene)^ 
which has a different dispersion relation than the one 
discussed here. Ref. [2^ presents a detailed analysis of 
these resonances for the case of a potential well, in the 
context of a two-dimensional massless Dirac equation. 
The way the plots of the transmission probability vs. in¬ 
cident energy in Ref. depend on the most important 
model parameters is a description of the Fano-type reso¬ 
nances which is complementary to ours. By calculating 
the Wigner delay they also predict the tempo¬ 

ral trapping that we are going to show below explicitly. 

In our present case a Fano-type resonance may occur if 
the energy of the scattered electron, after losing a number 
of quanta ntiv^ coincides with the bound state energy e of 
a Dirac particle in the time-independent potential, lying 
between —mc^ and mc^: 

Eq + nhv = e, n = —1, —2,... (30) 

In the following we focus on positive values of Vq, as it 
is depicted on Fig.[T] For L = I, one usually finds a single 
bound state of the static potential barrier, which simpli¬ 
fies the interpretation of the results. The energy eigen¬ 
values e corresponding to these states are in the range of 
[—mc^,mc^] , thus they are evanescent in regions 1 and 
3, where the potential is zero. Inspecting the positions of 
the resonances in Fig. [SJ we find that they appear when 

Eq + nhv ~ e, (31) 

with n being integer. Note that we cannot expect ex¬ 
act equality in Eq. (EH, since AC Stark shift modifies 
the energy levels. In the investigated parameter ranges. 




FIG. 8: Panel (a): {T) as a function of barrier height Vo for 
the indicated values of the oscillation amplitude D. Addi¬ 
tional parameters: Eq = = 0.2,1, = 1. Panel (b): 

Sideband transmission amplitudes tn as a function of barrier 
height Vo, for comparison with the D = 0.2 curve of panel (a). 
Note that the curves corresponding to n = —1 and n = —2 
have been rescaled (divided by a factor of 160). 



FIG. 9: Density plot of p{z,t), for L = 1, Vb = 3.97mc'^, 
Eq = l.lmc^ and D = 1.4. 


the relative difference of the left end right hand sides of 
Eq. (1311) was around 5% for low barriers (in the sense that 
Vo is not considerably larger than HVl), and it decreased 
for larger values of Vq ■ 

As one can expect, resonances corresponding to in¬ 
creasing magnitudes of n are weaker. According to the 














terminology that is often used - in spite of the fact that 
the potential oscillations are not quantized - we may say 
that the probability of higher order processes that involve 
2, 3, etc. excitation quanta (’’photons”) are considerably 
lower than that of ’’single photon” processes. However, 
when the amplitude of the oscillation increases, higher 
order resonances get more pronounced. 

Panel b) of Fig. |5]shows a clear example that whenever 
there is a resonance in (T ), the transmitted amplitude to 
corresponding to the incoming electron’s energy has a 
sharp minimum. At the same value of Vq, one of the 
coefficients (that belong to a;„ = Eq/H + nv, with 
n < 0) has a peak. We can clearly identify the peaks 
corresponding to n = —1 and —2 in Fig. Ob). 

As an example for a third order resonance, Fig.|9]shows 
p{z, t) for parameter values where t _3 has the highest 
magnitude amongst all transmission coefficients. As we 
can see, the fact that the incoming electron excites a 
localized bound state appears in this figure as an increase 
of the electron density. Note that this temporal trapping 
in region 2 is considerably more pronounced than in the 
cases seen in Figs. [Hand [T] the maximal value of p{z,t) 
inside the oscillating potential is five time larger than 
anywhere else. 

Figs. [5] and IHl demonstrate that the abrupt changes 
in the cycle averaged transmission and reflection proba¬ 
bilities are due to Fano-type resonances as described by 
Eq. (IMl) . Note that this kind of behavior can appear in 
principle for higher energies as well, but then the order 
of the process [the values of n in Eq. (IHTl) ] is so high, 
that it means only a practically invisible correction for 
the transmission probabilities. 

Finally, let us analyze to what extent the effects pre¬ 
sented so far are relevant for graphene. As we have men¬ 
tioned earlier, according to Ref. [s^ a band gap A as 
large as 0.26 eV can appear in epitaxially grown graphene 
on SiC substrate. This band gap plays the role of the en¬ 
ergy difference of 2mc^ that separates positive and nega¬ 
tive energy eigenstates of the massive Dirac equation we 
considered. Our findings are relevant when hfl, Vq and 
hv has the same order of magnitude as A. The charac¬ 
teristic frequency z//27r is in the THz regime {hv Ri A/2 


gives z/ = 25 THz), which is in the experimentally achiev¬ 
able range. That is, generation of wave packets, the ex¬ 
istence of alternating relativistic currents in the reflected 
region as well as the appearance of Fano-type resonances 
can be visible also in the case of the graphene. 


IV. SUMMARY 

We presented a theoretical work on the relativistic 
scattering of massive Dirac particles on a time-periodic 
rectangular potential barrier, using Floquet theory. We 
used the dependence of the cycle-averaged transmission 
probability on the barrier height to describe the quasis¬ 
tationary behavior of the system, in the case of a weakly 
relativistic, relativistic and ultra-relativistic incident par¬ 
ticle. In each case, the system shows Klein tunneling 
and Fano-type resonances. We explored the details of 
the transport with the help of space and time dependent 
currents and densities, which show explicitly that the os¬ 
cillating barrier generates wavepackets from the incident 
plane wave. We explained in detail the Fano-type res¬ 
onances with the interplay of the sideband states gen¬ 
erated by the oscillating potential and the bound states 
of the barrier, and showed the corresponding temporal 
population trapping in the barrier region. Finally we 
discussed the relevance of our results to graphene with 
an induced bandgap. 
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